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ABSTRACT 

0> ' We calculate how accurately parameters of the short-period binaries 

^ ' (10"'* Yiz ^ / ^ 10~^ Hz) will be determined from the gravitational waves by 

^ LISA. In our analysis the chirp signal / is newly included as a fitting parameters 

f — ' and dependence on observational period or wave frequency is studied in detail. 

oo ■ 

^ ' Subject headings: gravitational waves - gravitation - binaries 

o ■ 
o 

eg ■ The Laser Interferometer space Antenna (LISA), a joint project of NASA and ESA 

(European Space Agency) would establish gravitational wave astronomy at low frequency 



Implications for gravitational wave astronomy are also discussed quantitatively. 



1. Introduction 



band (10~^ H2; ^ / ^ 10~^ H2;). It would bring us essentially new information of the 
Universe (Bender et al. 1998). For example gravitational waves from merging super massive 
black holes (SMBHs) would be detected with significant signal-to-noise ratio [SNR) > 10^, 
though event rate of such merging is highly unknown {e.g. Haehnelt 1994, Vecchio 1997). 
Galactic binaries are promising sources of LISA (Mironowskii 1965, Evans, Iben & Smarr 
1987, Hils, Bender & Webbink 1991, Webbink & Han 1998). Gravitational waves from some 
known binaries {e.g. X-ray binary 4U1820-30) would be detected with SNR > 5 by one 
year integration (Bender et al. 1998). In addition more than thousands of close white dwarf 
binaries (CWDBs) are expected to exist in LISA band (for recent studies see Yungelson et 
al 2001, Nelemans et al. 2001, Napiwotzki et al. 2002). Our target in this article is these 
Galactic binaries. We examine how accurately information of binaries can be extracted from 
gravitational waves observed by LISA. 

Cutler (1998) studied the estimation errors for binary parameters with special attention 
to angular variables such as direction and orientation of binaries (see also Peterseim et al. 
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1997; Vecchio & Cutler 1998; Cutler & Vecchio 1998; Hughes 2002; Moore & Hellings 2002). 
He used approximation that emitted gravitational waves would be monochromatic, namely 
neglected the effects of the chirp signal /. But the wave frequency or the chirp signal are 
fundamental quantities for gravitational wave astronomy. From the measured chirp signal 
/ wc can obtain the so called chirp mass Mc = Mf^^M2^^{Mi + M2)~^^^ (Mi,M2: masses 
of two stars) for a binary whose orbital evolution is determined by gravitational radiation 
reaction. Furthermore the distance to the binary could be estimated from the chirp signal 
/ and the amplitude of the wave signal (Schutz 1986). The frequency / itself also contains 
important information. One of the authors (Seto 2001) pointed out that signature of the 
periastron advance could be detected in gravitational waves from an eccentric binary by 
measuring its wave frequencies preciously. If this method works well, we can estimate the 
total mass Mtotai = + ^2 of the binary beside the chirp mass, and consequently each 
mass of the binary is obtained separately. 

Estimation errors for fitting parameters correlate complicatcdiy to each other and de- 
pend largely on observational situations. For example longer observational periods would 
improve not only signal to noise ratio but also resolution of the frequency space. Note that 
the latter is crucial for reducing Galactic binary confusion noise, as the number of resolved 
binaries increases with decrease of the frequency bin oc T^l (see Seto 2002 for details). At 
LISA age it would become an interesting observational challenge to optically identify the 
binaries whose gravitational waves are detected by LISA. We also discuss impacts of these 
observational efforts on estimation of binary parameters. 

The magnitude of the chirp signal / has strong dependence on the wave frequency / as 
/ oc mI^^ f^^l^. At lower frequencies the estimation error A/ would be larger than the signal 
/ itself. Then it would be better to remove / from fitting parameters and simply put f — 
from the beginning. By evaluating the estimation error A/ we can quantitatively discuss 
these prescriptions for signal analysis. 

In this article we only study parameter estimation errors on the assumption that (i) 
signal has been detected with high SNR, and (ii) noises are Gaussian distributed. We evaluate 
the estimation errors using the Fisher information matrix for maximum likelihood method. 
For low SNR our simple analysis would not be valid, as the probability distribution function of 
the fitting parameters could become highly complicated {e.g. multimodal) (Balasubramanian 
& Dhurandhar 1998). Thus our result should be regarded as a lower bound of the estimation 
errors. Non-Gaussian nature of noises would farther increases errors. This is a very important 
problem, but detailed quantitative analysis would be a formidable task. 

This article is organized as follows. In §2 wc briefly discuss the gravitational waveforms 
of chirping binaries and the data stream obtained by LISA. Then we mention the parame- 
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ter estimation based on the Matched filtering analysis. In §3 we numerically evaluate the 
parameter estimation errors, and discuss its dependence on observational period or wave 
frequency in detail. §4 is devoted to summary and discussions. 



2. Gravitational Waveforms and Parameter Extraction 

2.1. Gravitational Wave Measurement with LISA 

LISA consists of three spacecrafts forming an equilateral triangle, and orbits around 
the Sun, trailing 20° behind the Earth. The sides of the triangle arc L = 5 x 10^ km 
in length, and the plane of triangle is inclined at 60° with respect to the ecliptic. The 
triangle rotates annually. The gravitational wave signal is reconstructed from the three data 
streams that effectively correspond to three time- varying armlength data, {SLi, SL2, SL3) for 
gravitational waves with X < L. We basically analyze two data streams given by si{t) — 
{SLi{t) - SL2{t))/L and sn{t) = {SLi{t) + SL2{t) - 25L^{t))/ ^/2,L. These data can be 
regarded as the response of two 90°-intcrfcromcters rotated by 45° to one another (Cutler 
1998). The data siji{t) contain both gravitational waves signal hiji(t) to be fitted by 
matched filtering and the noise niji{t). The latter is constituted by the detectors noise and 
the binary confusion noise. As in Cutler (1998) we assume that the noises are stationary, 
Gaussian and uncorrelated with each other. 

The gravitational wave signals hijiit) from a binary are written as 

where F^\^ [t) are the pattern functions which depend on the source's angular position of the 
binary (ds, 4>s) , its orientation {9l, and detector's configuration. The angular variables 
with bars are defined in a fixed barycenter frame of the solar system, and the direction and 
orientation of the binary are assumed to be constant during the observation in this frame. 
The quantities /i+^x(^) are the two polarization modes of gravitational radiation from the 
binary. We can estimate both {9s, ^5) and {9^, (pi) from the time profiles of the two signals 
due to LISA'S annual rotation and revolution. Further discussion and details about the 
pattern functions arc seen in Cutler (1998). We basically use his formulation but newly 
include the effects of the chirp signal /. 
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2.2. Gravitational Waveforms 

We study short-period (10^'' Hz < / < 10^^ Hz) binaries such as the close white dwarf 
binaries (CWDBs) or the neutron star binaries (NBs) in our Galaxy (Hils, Bender & Web- 
bink 1991). We only discuss binaries with circular orbits. This is an excellent approximation 
for CWDBs, as their orbits are circularized by strong tidal interaction in their earlier evo- 
lutionary stages. At LISA band some NBs or neutron star-white dwarf binaries would have 
non- negligible eccentricities e ~ 0.1 induced at their supernova explosions (see e.g. Brown, 
Lee, Zwart, & Bethe 2001). But extension to these eccentric binaries is straightforward. 



The chirping gravitational waveform is given by the quadrupole approximation (Peters 1964) 



as 



h+it) 

Kit) 



A COS 
-2A sin 



X 



1 • 



2Tl[f + -ft]t + 4>D{t)+4>0 



1+ (l-u^' 
X ( L ■ n ) , 



(2) 



where L (given by Ol, 4>l) is the unit vector in the direction of the binary's orbital angular 
momentum, fi (given by 9s, 4's) is the unit vector toward the binary and 0o is an integral 
constant. We regard the frequency / and its time variation / as constants in the above 
equations. A purely monochromatic waveform has / = 0. This is the case studied by Cutler 
(1998). When gravitational radiation reaction dominates evolution of the binary as in the 
case of CWDBs or NBs, the chirp signal / is given as / = (967rS/V5)/ii/3Mc/^ with the chirp 
mass Mc- The perturbative expansion for the intrinsic evolution of the wave frequency in 
equation (2) are valid under the condition fTobs <^ /, where Tabs is the observational period. 
This condition is expressed as 
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The amplitude A in equation (2) is given in terms of the wave frequency /, chirp signal / 
and the distance D as 

A = — L_ (4) 

Thus we could determine the distance if we could measure three observables /, / and A 
(Schutz 1986). This is an important aspect of gravitational wave astronomy. 

The term 0_d(^) in equation (2) is caused by revolution of LISA around the Sun and 
called the Doppler phase. Its explicit form is given by 



(t)D{t) = 27r/i?sin^5COS [0(i) - ^5] , 



(5) 
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where R—1 AU and (f){t) — 27rt/T {T — lyr) is the direction of LISA in the fixed barycenter 
frame. 



2.3. Parameter Extraction 

Let us briefiy discuss the matched filtering analysis and the parameter estimation errors 
(Finn 1992; Cutler & Flanagan 1994). We assume that the signal ha{t) is characterized by 
some unknown parameters 7^ (eight parameters in the present case: {A, f, f, 0o, ^s, 4's, ^l, 
In the matched filtering analysis the variance- covariance matrix of the parameter estimation 
error A'ji is given by inverse of the Fisher information matrix Tij as (AjiAjj) = (F"^)^^.. 
For a quasi- monochromatic binary {fTobs ^ /) the noise spectrum Sn{f) is nearly con- 
stant in the frequency region swept by the binary and the Fisher matrix simply becomes 
(Cutler 1998) 



The error boxes for the angular parameters L and fi become ellipses in the celestial 
sphere due to the correlation of two parameters 9 and 0. In this article we represent the 
estimation errors for direction and orientation of binaries in the form defined in Cutler (1998) 
as follows 

An = 27rV(A/x')(A02) - (A//A(/.)2, (7) 

where we have defined = cos^. In the same manner the signal to noise ratio (SNR) is 
given by 

(SNR) ' = Yl I °" K{t)K{t). (8) 

From equations (6) and (8) it is apparent that the expressions for the estimation errors 
(A7jA7j) do not depend on the noise spectrum Sn{f) when they are normalized by the signal- 
to- noise ratio (Cutler 1998). In this article we extensively use this normalization method. 

For example the parameter estimation errors for a simple wave form h{t) = Asin[27r(/ + 
ft/2)t + 0o] with four fitting parameters {A,f,f,(l)o) is easily evaluated as in Seto (2002). 
The explicit forms for AA, Af and A/ are given by 

A SNR ^ [ 10 ; ' ^ ^ 

_ l^i- 22 r^V'T- (10) 
" TT SNR-^-^\ 10 J 
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This simple analysis does not include information of the angular parameters, but would be 
helpful to understand more detailed numerical analysis in the following section. 



3. Results 

3.1. General Behavior 

We have numerically evaluated the uncertainties of the estimated parameters for various 
quasi-monochromatic binaries. In this subsection we show results for a typical example 
with a fixed set of angular parameters at cos ^5 = 0.3,05 — 5.0, cos = —0.2 and (pL = 
4.0. In Table 1 wc show LISA's measurement accuracy for parameters {A, f, f ,Qsi^L) 
at frequencies / = 10^'', 10^'^ and lO^^Hz. We present our results for two observational 
periods, Tq^^ = 1 and lOyr. All results are normalized by SNR = 10 after integration over 
each observational period Tgbs- These results simply scale as {SNR/10)~^ for errors AA, Af 
and Af and {SNR/10)~'^ for error ellipses AO,s,l- The first row for each observational period 
Tabs represents the case when all the eight parameters {A, f, f, (po, 9s, 05, Ol, 4>l) are included 
in the matched filtering analysis. The second row corresponds to the case without the chirp 
signal /. Results for Tobs — lyr are obtained under the same condition with Table 1 (case A) 
in Cutler (1998). Our numerical values completely coincide with his ones. The third row is 
the case when direction of the binary (^5, cps) is given exactly by other method {e.g. optical 
identification of the binary) and removed from the fitting parameters. 

Figs 1 and 2 show clearly dependence of the orbital parameters on wave frequency and 
observational period Tobs- In Fig-1 we show LISA's measurement accuracy as a function of 
observational period Tobs at given frequencies / = 10^'^, 10^'^ and 10^^ Hz. All results are 
normalized by SNR = 10 at integration period Tobs = 1 yi'- The solid lines are results for 
fitting all the eight parameters. The dotted lines represent the case when the angular position 
{^s,4>s) are removed from the fitting parameters (see also Hughes 2002). The dashed fines 
are results with fitting only the angular position {9s, 4>s)- For observational period Tgbs ^ 2 
yr, the difference between the sofid and the dotted fines is very small especially for AA, Af 
and Af irrespective of the frequency. Thus optical determination of the source direction fls 
would only slightly reduce the estimation errors of other parameters. Asymptotic behaviors 
of errors AA, Af and Af are given by 



(12) 
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Af = 0.22 (13) 

A/ = 0-43 (^)"t-, (14) 

(see also Table. 1), and the asymptotic time dependence is given as AA oc T^^y^, Af oc T'^J"^ ^ 
and A/ oc T'^J"^ , since we have S'A'"i? oc tH"^ from equation (8). We have also AQ:s,l oc 
SNR-^ oc T'oftg for directions and orientations of binaries. Numerical results for Af and 
Af in equations (13) and (14) are almost identical to analytical ones in equations (10) and 
(11). Hence we can expect that the information of the angular parameters (the direction and 
orientation) does not affect the accuracy of the estimation for the frequency / and the chirp 
signal / for observational period T^bs ^ 2 yr. But for the amplitude A, the above result (in 
[12]) is two times as large as equation (9). This is due to the fact that the amplitude A is 
tied with the angular parameters (the inclination L • n) in equation (2) and the estimation 
error AA strongly depends on the error in L • n (see the detailed discussion in the next 
subsection) . 

Fig.2 is same as Fig.l but given as a function of frequency / with fixed observational 
period Tobs = 1 and 10 yr. We find that for Tobs = 10 yr the errors AA, Af and Af do not 
depend on the frequency, as given in equations (12)-(14). But the angular resolution AQs 
depends on the frequency and is given by 

AQs = 4.8xl0-^f— ^VV^V'sr, (15) 

for higher frequencies / > 10~^ Hz, and nearly constant for lower frequency / < 10~^ Hz. 

Now we discuss correlation between Fisher matrix elements for T > 2yr. As expected 
from Figs 1 and 2, the source direction has almost no correlation with other parameters at 
higher frequencies / > 10~^Hz ^ and weakly correlates with source orientation at lower 
frequencies / < 10~^Hz. This dependence seems reasonable considering the information 
used to determine the source direction Qs, namely, the Doppler phase at / > 10~'^Hz and 
amplitude modulation at / < 10^'^Hz. The orientation is strongly related to the latter. 
As a result angular variables Qs arid Ql correlate at lower frequencies. The detected wave 
form ha is modulated by annual revolution and rotation of the detectors, and its phase 



^In other words, Fisher information matrix Tij in equation (6) becomes diagonalized to the two parts; 
r ~ ( -n* ) ' where Ta represents a 6 x 6 matrix which corresponds to {A, /, /, cpo, 9l, (j)^) components 



, Ts, 

and Vs represent a 2 x 2 matrix which corresponds to {ds, (j^s) components. 
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has oscillating part due to the modulation that is more prominent at higher frequencies. 
Amplitudes of the derivatives dh/dQg also show this oscillation. In the Fisher matrix their 
cross terms with other derivatives are significantly canceled with long time integration. As 
a result the estimation errors for the angular direction fig show very weak correlation with 
other errors. 

In the same manner wc can also understand the asymptotic frequency dependence of the 
angular resolution O^. In the derivatives dh^/dji for the Fisher matrix elements in equation 
(6) the frequency / appears only through the Doppler phase (pnit) in equation (5). At lower 
frequencies / < terms from the Doppler phase (oc /^) is smaller than the terms 

from the amplitude modulation (oc Z*^), and consequently the estimation error AQs does 
not depend on the frequency /. At higher frequencies / > 10~^Hz the Doppler phase term 
becomes dominant and the Fisher matrix are diagonalized. Thus the error AQ5 depends on 
the frequency / as AQ5 oc f'^ (Cutler & Vecchio 1998, Moore & HeUings 2002). 

Finally we analytically investigate how the estimation errors improve with the observa- 
tional period Tof,s- We reanalyze the simple toy model h{t) = y4sin[27r(/ + ft/2)t + 0o] (see 
the sentences before equation (9)) adding the information of the source direction {9s, 4>s) by 
the Doppler phase 0d(^)- This waveform is given as 

h{t) = A sin[ 2n{f + ft/2)t + Mt) + 0o], (16) 

with six fitting parameters 7^ = {A, f, f , (f)o,9s, (ps)- We evaluate the magnitude of the 
variance-covariance matrix {A'fiA'-fj) by its determinant as det(A7jA7j) = (detF)"^. After 
some algebra the time dependence of detF is given analytically from equations (6) and (16) 
as, 

.12 6 \ 90 



det F oc 1 2 1 J , (17) 

V TT^ (r„,,/lyr)V V 7rMWlyr)V ^ ' 

where we assume that the observational period Toi,s{^ f~^) is a integer in units of year. 

In the above expression the quantity (detF)~^ formally diverges at Tabs = v'^/vr = 0.98 

[yr] which is very close to 1, thus the estimation errors have large value at Tabs = 1 yr- 

Fig. 3 shows the inverse of the determinant of the Fisher matrix in the form Tj£,(det F)"^ as 

a function of the observational period T^bs- From Fig.3, Tj^^^(det F)^^ rapidly converges at 

Tobs ~ 2 yr. In this case we also obtain the estimation errors as follows; 

AA 1 

-A - SNR' ^^^^ 



_ 4^3^^ l-(45/8) x^{l+x-) 

~ TT SNR ]l (l-6x2)(l-90x4) ' ^ 

■ _ 6V5 t;,1 1 

- TT SNR ^/TTgo^' ^^"^ 
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^^s = — = ^ X , (21) 



_4 /cos^y' / / fSNR\-^ 

V 0.3 ) Vio-2Hzy V 10 / 



^ 4.3x10-^:::^ [—^\ sr, (22) 



where x = 7r~-'^(T'o6s/lyr)~''^ and equation (22) is valid for To^g > 2 yr and is similar to 
equation (15). The above results for A/, A/ and AQ5 are about 40% smaller than the 
the results in Fig.l at Tobs — 1 yr, because we do not include the information of the source 
orientation. A/ and A/ are asymptotically same as equations (10) and (11) for Tabs ^ 2 yr. 



3.2. Statistical Analysis 

So far we have studied a specific set of the angular parameters ft and L. In this subsection 
we present statistical results for their various combinations at the asymptotic region Tobs ^ 
2yr. 

We have made 100 realizations of h and L that are distributed randomly on celestial 
spheres. Then we calculate the estimation errors A/ and A/ for each binary normalized 
by SNR — 10 with Tgbs — lOyr. We find that theses errors depend very weakly (less than 
10% scatter) on the directions n and L in contrast to the results for T^bs — lyr (scattering 
typically factor 2=*=^ around the average). We calculate their mean values and obtain results 
given in the following forms 

A/ = 0.22 (^)"'t-, (23) 

A/ = 0-«(^)"r,iJ- (24) 

Note that these results do not depend on the frequency / and would be useful for quantitative 
analysis of quasi-monochromatic binaries with Tgbs ^ 2yr. As expected from the previous 
subsection the simple analytical estimations given in equations (10) and (11) are very close 
to equations (23) and (24). 

We have also studied the estimation errors AQ5, AJl^ and AA/A. Distributions of the 
latter two have very large scatters. This is because (i) the Fisher matrix elements relating 
to the orbital orientation L become singular at the highly symmetric face-on configuration 
n — ±L and (ii) the estimation of the amplitude A is closely related to the inclination. For 
majority of realizations with \n • L\ ^ 0.8 we have typically AA/A ~ 0.2(5'iVi?/10)^^. To 
determine the source direction n we can use the information of the Doppler phase in addition 
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to the annual amplitude modulation caused by LISA's rotation. Thus the error AQ5 does 
not show such a bad behavior. For the above 100 realization at / = O.OlHz with SNR = 10 
and Tf^bs = 10 yr the errors AQs are distributed as 1.3 x 10~^sr < AQs ^ 3.7 x lO^'^sr 
with the mean value AQs = 7.1 x 10^'^sr. Therefore the following relation roughly gives the 
estimation error for the source direction with Tobs ^ 2yr 

Ans ~ 7.1 X 10-^ ( ( TT^V'sr, (25) 



at / > 2 x 10 ^Hz where the Doppler phase becomes more important than the annual 
amphtude due to the rotation of LISA. 

Now let us discuss issues concerned with the chirp signal / based on equation (24). We 
assume that the chirp signal is dominated by gravitational radiation reaction and introduce 
a parameter R = A/// that represents relative accuracy of the measured signal /. For a 
given threshold R we can solve the corresponding frequency fa as 

for observational period Tghg > 2yr. The frequency Jr^i can be regarded as the critical 
frequency for treatment of the chirp signal /. At higher frequencies / > fR=i the parameter 
/ should be included in the matched filtering, but the simple prescription / = would be 
better at / < /i?=i since the expected signal would be completely buried in error. If the 
chirp signal is measured with accuracy i?(-C 1), the chirp mass can be estimated with relative 
accuracy AMc/Mc ~ 3R/5 from relation / oc Mc^^. The chirp signal / is also essential to 
determine the distance D to the binary. From the simple relation A — 5//967r^/^£), the 
estimation error for distance D is roughly evaluated as 

AD 
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1.4xlO-3Hzy \IMqJ VlOyr 

where we have used the typical error for the amplitude estimation AA/A ~ 0.2(5'A^i?/10)^^. 
For compact binaries with chirp mass Mc ~ IM© such as NBs or CWDBs the chirp signal 
/ is the dominant source of the error at lower frequencies / < 10~^Hz and observational 
period Tofes ~ lOyr. The amplitude A becomes dominant one at / > 10~^Hz. 
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Finally, we calculate the estimation error for the three dimensional position of the 
Galactic binaries. The signal to noise ratio is calculated from equation (8) as SNR ~ 
380 (//5xlO-3Hz)2/3 (D/10kpc)-i (T,5,/10yr)i/2 {MjlM^f/^ (v^/8 x 10-2iHz'^/2)-i/2^ 

The noise spectrum Sn{f) (in units of Hz^^) is nearly constant for 3 x 10^'^ Hz < / < 10^^ 
Hz (see Fig.5 in Cutler 1998). Hence one could determine both the angular position and the 
distance to the binary with the accuracy of 



-5/3 / r, \ / ^ \ -V2 



5xlO-3Hzy VlOkpcy 



lOyr, 

-5/3 / fn- \ 1/2 



"<litj Uxl0-Hz-V2) ™' ^30) 



5xl0-3(_ .^„V^7^^ (^V^' 



AD 

H^^J U X 10-21HZ-V2 J • ^^^^ 



4. Conclusion 

We have calculated LISA's measurement accuracy for short-period binaries (10~^ Ylz < 
f < 10~^ H^) in our Galaxy, including the effects of chirp signal / and dependence on 
observational period Tgbs- We find that the measurement accuracy rapidly improves for 
Tobs ^ 2yr comparing with To^s ~ lyr- This might be an important element for discussing 
operation period of LISA. 

At observational period Tobs ^ 2 the errors for quantities / and / is independent on 
the information of the positions and orientations of binaries in contrast to AA and AQs,l- 
It is also found that the estimation errors AA, Af and Af are almost independent on the 
frequency /. The fitting formulae of the estimation errors for the source parameters (such as 
frequency /, chirp signal /, amplitude A, angular position Qs and distance D) are given as 
functions of / and Tobs- We expect these would be powerful tools for quantitatively studying 
possibility gravitational wave astronomy. 

The relative motion between the Sun and a binary could affect the frequency / by the 
Doppler factor, (I + vr) where vr is the relative radial velocity. The tangential velocity vt or 
accelerating motion v would affect the time variation of the frequency, {f / f)gai ~ v^/cD or 
v/c v^ot/cRg respectively where Vrot is the rotation velocity (~ 200km/s) of the Galaxy, 
Rg is the Galactic radius ~ lOkpc and D is the distance to the binary (see e.g. Damour 
& Taylor 1991). For most Galactic binaries in LISA band this should be negligibly small 
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compared with {f / f)GW due to gravitational radiation reaction. This condition is expressed 
as / > 3 X lO-^(Mc/lM0)-5/8(t;rot/2OOkm)3/4(L)/lOkpc)-^/s Hz. For binaries very close to 
VLS D < 100 pc these effects would be important. But note that we might measure the proper 
motion of such binaries by optical observation. 
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manuscript, and Takeshi Chiba, Hideyuki Tagoshi and Hirotaka Takahashi for useful com- 
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Table 1: LISA's measurement accuracy for binaries with angular parameters (cos ^5 = 
0.3,05 = 5.0, cos = -0.2, 0i = 4.0). Results are normalized by SNR = 10. Errors 
scale as (5iVi?/10)-^ for AA, A/ and A/, and {SNRjVS)-'^ for A^s^l- The second lines 
in each observational period Tobs represent the case with removing the chirp signal / from 
fitting parameters and the third lines with removing the direction of the source (^5, ^5). 
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Fig. 1. — LISA'S measurement accuracy for the binaries as a function of the observational 
period Tobs with angular parameters (cos^^ = 0.3,05 = 5.0,cos^^l = — 0.2,0^ = 4.0). The 
solid lines correspond to A/, A/, A A, AQl and Afls from top to bottom. The dotted lines 
represent the case source positions {9s, 4>s) are removed from the fitting parameters, and 
the dashed lines represent the case the only source positions are included in the fitting 
parameters. The accuracies are normahzed by SNR — 10 at lyr observation. 
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Fig. 2. — Same as Fig.l, but as a function of the frequency / with the observational period 
Tabs = 1, 10 yr. The accuracies are normahzed by SNR = 10. 
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Fig. 3. — The inverse of the determinant of the Fisher matrix for the simple waveform given 
in equation (16). The filled circles represent with integer Tgi^g in units of year. We normalize 
the overall scale by unity at Tobs — +oo. 



